Probing small-scale intermittency with a fluctuation theorem 
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We characterize statistical properties of the flow field in developed turbulence using concepts from 
stochastic thermodynamics. On the basis of data from a free air-jet experiment, we demonstrate 
how the dynamic fluctuations induced by small-scale intermittency generate analogs of entropy- 
consuming trajectories with sufficient weight to make fluctuation theorems observable at the macro- 
scopic scale. We propose an integral fluctuation theorem for the entropy production associated with 
the energy transfer along the eddy-hierarchy and demonstrate its extreme sensitivity to the accurate 
description of the tails of the velocity distributions. 
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All processes in nature are bound to produce entropy. 
This central dogma of macroscopic thermodynamics got 
substantially qualified in the preceding decade by new 
insights into the properties of small, strongly fluctuating 
systems. If entropy consuming trajectories occur with 
appreciable probability, thermodynamic inequalities may 
be considerably tightened to assume the form of equalities 
[3, 0] • The emerging field of stochastic thermodynamics 
(for recent reviews see 0, 0] ) focuses on the full proba- 
bility distributions of thermodynamic variables like heat, 
work, and entropy and establishes thermodynamic rela- 
tions for individual fluctuating histories of the systems 
under consideration. Most prominent among these re- 
lations are the so-called fluctuation theorems quantifying 
the relative frequency of entropy-consuming as compared 
to entropy-producing trajectories. Applications of these 
developments concern free-energy estimates of biopoly- 
mers p, [1| , the efficiency of nano- machines 0, Hj , and the 
thermodynamic cost of information processing tlj fioj , to 
name a few. 

On the experimental side, most investigations have 
been done with nanoscopic setups like in single-molecule 
manipulations IT 13j, colloidal particle dynamics [l^ - 
flH | or harmonic oscillators [l6j]. For these systems, typ- 
ical free-energy differences are of order k^T and the 
ubiquity of thermal fluctuations ensures the broad dis- 
tributions of work, heat, and entropy which are indis- 
pensable for the application of fluctuation theorems. In- 
creasing the size of these systems to macroscopic orders, 
the importance of thermal fluctuations fades, entropy- 
consuming trajectories become exceedingly rare, and 
the fluctuation theorems degenerate to the inequalities 
known from traditional thermodynamics. Besides some 
investigations in granular media [l7l - fl9j ]. rather few ex- 
amples of macroscopic systems have been identified which 
are amenable to an analysis using fluctuation theorems. 

Turbulent flow of liquids and gases is a fascinating phe- 
nomenon with many different facets that has been cap- 
tivating scientists for centuries. Despite its broad range 
of technical relevance including turbulent drag [2^ , tur- 
bulent mixing plj , atmospheric turbulence with implica- 
tions for climatic models 22( and the prospects of wind 




FIG. I: Distribution p(u,r) of velocity increments u at vari- 
ous scales r (circles) in the turbulent flow of a free jet exper- 
iment [24|. The velocity increments u are given in units of 
the standard deviation — 0.54 m/s at infinite scales. Also 
shown is the instantaneous stationary distribution p st (u,r) 
defined in ((Sj) (full lines) and Gaussian fits to the experimen- 
tal data (dashed lines). Both the deviation from the Gaussian 
approximation and from the stationary distribution increases 
when approaching smaller scales. For the sake of clarity, the 
distributions for various scales are vertically shifted by 10 3 . 



energy [23|, several aspects of turbulent flows are still 
not fully understood. In particular, the intricate pattern 
of small-scale flow in developed turbulence with its in- 
termittent change between laminar periods and violent 
bursts of activity have eluded a thorough theoretical un- 
derstanding so far. 

In the present letter, we show that the fluctuating flow 
field of developed turbulence represents a proper test sys- 
tem for stochastic thermodynamics. The dynamic fluc- 
tuations of turbulence show up at a macroscopic scale 
and, at the same time, are strong enough to generate 
"non-mainstream" trajectories with sufficient frequency 
to observe fluctuation theorems in action. Using data 
from a free air jet experiment, we elucidate the nature 
of the entropy-consuming trajectories and demonstrate 
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the convergence of an integrated fluctuation theorem for 
data sets of rather small size. We further discuss how to 
use the fluctuation theorem for the statistical description 
of the flow field. 

Previous applications of fluctuation theorems to turbu- 
lence have been restricted to versions of the steady-state 
fluctuation theorem [l], [25[ . Experimentally, fluctuations 
of heat flux j2(| and pressure [27| have been studied; nu- 
merical investigations concerned the fluctuations of in- 
jected 
tion 13: 
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29j and hydrodynamic entropy produc- 
. In contrast, we propose an integral fluctuation 
theorem for the total dynamic entropy production char- 
acterizing the non- stationary energy transport in the in- 
ertial range of developed turbulence. 

In a standard setup, isotropic turbulence is generated 
by injecting energy into the flow by an external force 
field at a large, so-called integral scale L [3l], [32jj. By 
repeated break-up of eddies, a self-similar eddy hierar- 
chy forms which is characteristic for developed turbu- 
lence [33^ . On average, energy is transferred along the 
cascade to smaller and smaller scales until, due to molec- 
ular friction, it is dissipated in the viscous range. The 
Taylor scale A marks the length scale above which the 
influence of dissipation is still negligible. 

A suitable way to characterize the stationary, homoge- 
neous, and isotropic flow field v(x, t) in the inertial range 
between L and A is via the probability densit y fu nction 
p(u, r) of longitudinal velocity increments [3lL l32| 



u(r) 



(v(x + er, t) — v(x, t)) 



(1) 



Here, r denotes the scale at which the velocity difference 
u is evaluated, e is a unit vector and due to the average 
symmetries of the turbulent flow, the statistical proper- 
ties of u only depend on r. Fig. Q] shows histograms of 
this distribution using data obtained in a turbulent air jet 
experiment In this setup, L = 6.7 cm, A = 6.6 mm, 
and the nozzle-based Reynolds number is about 2.7- 10 4 . 
The flow velocity v(t) is measured a distance of 125 nozzle 
diameters away from the nozzle and then converted to a 
flow field v(x) by use of the Taylor hypothesis. Chopping 
v(x) into non-overlapping intervals, N = 5 • 10 trajecto- 
ries u(r) are obtained from which the shown histograms 
are compiled. As Fig. [1] clearly shows, p(u, r) exhibits a 
Gaussian form for scales r ~ L and develops pronounced 
non-Gaussian tails towards scales r ~ A. This effect is 
commonly referred to as small-scale intermittency, as in- 
termittent bursts in v(x) cause the boosted occurrence 
of large values of u on small scales [IH • 

An inventive approach to characterize the properties of 
the distribution p(u, r) in the inertial range is to interpret 
u{r) as realizations of a Markov process on the eddy hier- 
archy with the scale r playing the role of time [35j . The 
evolution of p(u, r) is then described by a Master equa- 
tion with initial condition at r = L, for which a Kramers- 
Moyal expansion (36j may be performed. For a variety 
of experimental situations, the Markovian character of 



u(r) was verified, and the coefficients of the corre- 
sponding Kramers-Moyal expansion were determined on 
the basis of experimental data 37-4(|. Moreover, in the 
limit of large Reynolds number, it is possible to systemat- 
ically derive the Master equation governing the evolution 
of p(u, r) from the underlying Navier-Stokes equations of 
the fluid flow and to recursively calculate the coefficients 
[il], HI]. In either way, one finds that drift and dif- 
fusion coefficients, and D^ 1 respectively, have well- 
defined, non-zero limits, whereas all higher coefficients in 
the Kramers-Moyal expansion vanish asymptotically. We 
are thus left with a Fokker-Planck equation (FPE) of the 
form 



-d r p(u,r\u L ,L) = [- d u D (1) (u,r) 



d 2 u D^(u,r)]p(u,r\u L ,L) 



(2) 



ruling the statistics of velocity increments on the eddy- 
hierarchy of developed turbulence. The minus sign on the 
l.h.s. of the FPE indicates that the evolution proceeds 
from large to small scales. 

The drift and diffusion coefficients, D^> and D^ 2 \ typ- 
ically depend on r and u; for the data shown in Fig. Q] 
one obtains, e.g., [24[ 

D^(u, r) = -a r - 6 - ai r~™ 7 u + a 2 u 2 - a 3 r a V (3) 

£>(*> ( u> r ) = 6 or O-25 _ bir 0.2 u + fc 2r -0.73 u 2 (4) 



with 

a = 0.0015, ai = 0.61, a 2 = 0.0096, a 3 = 0.0023 , 
6 = 0.033, 61 = 0.009, b 2 = 0.043 . 

The stochastic dynamics defined by |(5J) therefore exhibits 
characteristics of a driven non- equilibrium system. This 
is apparent also from the difference between p(u, r) and 
the instantaneous stationary distribution of the FPE @ 
for fixed scale r given by 

P st (", r) = , Z(r) = J e-*M du (5) 

with the stochastic potential 

V (u,r) = In (u,r) - j ^gr^d.' . (6) 

—00 

Examples of p st (u,r) have been included into Fig. [1] 

In the spirit of stochastic thermodynamics Q , we now 
associate with every individual trajectory u(r) a total 
entropy production 



AStot [«(•)] = _ / d r u(r)d u (p(u(r),r) dr 

i (7) 
p(ul,L) 
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FIG. 2: Empirical average (exp — A5tot)jv defined in ((9} for 
the experimental data of Fig. [1] as a function of the sample 
size N. According to the fluctuation theorem p]>. the aver- 
age has to converge to the horizontal line. The inset depicts 
the corresponding distribution of the total entropy production 
AStot as defined by (J7J). 



In the usual thermodynamic setting, the first term on 
the r.h.s. of ([7]) would describe the heat exchange with 
the reservoir, whereas the second one gives the entropy 
change of the system itself. The total entropy production 
([7]) fulfills the integral fluctuation theorem 43 1 



/«(•) 



1 



(8) 



where the average is over the different realizations of u(r). 

A reliable estimate of the exponential average in (JSJ) on 
the basis of a finite sample set is possible only if trajecto- 
ries with AStot [ u {')\ < occur with sufficient frequency. 
In the present setting, such trajectories correspond to re- 
alizations of the flow field that, contrary to the average 
tendency, transfer energy uphill the eddy-hierarchy We 
have used subsets of size N of the realizations for u(r) 
underlying Fig.[T]togcthcr with their entropy productions 
determined by (UJ and calculated the empirical average 
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1 N 
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_ A <JW 



(9) 



corresponding to ([5]). The results shown in Fig.[2]demon- 
strate that convergence to the asymptotic value is rather 
fast. This is corroborated by the appreciable weight of 
trajectories with negative entropy production in the dis- 
tribution p( AStot) shown in the inset. The macroscopic 
fluctuating flow fields of developed turbulence therefore 
share important features with the thermodynamic vari- 
ables of nanoscopic non-equilibrium systems under the 
influence of thermal noise. In particular, in both cases the 
respective probability distributions arc sufficiently broad 




FIG. 3: Typical form of measured velocity increments u(r) 
(full lines) realizing a very small (top) and a very large (bot- 
tom) entropy production AStot defined by ([7J. The dashed 
lines show the average part of v(x) neighboring «(•). The 
Taylor scale A and integral scale L are indicated by vertical 
lines. 



to allow an application of the concepts of stochastic ther- 
modynamics. 

The convergence of the empirical average to the 
theoretical value 1 given by © also indicates that the 
drift and diffusion coefficients ([3]), (jlj), estimated on the 
basis of the experimental data, describe the stochastic 
properties of the process u(r) rather well. Conversely, 
by monitoring @ during the numerical estimation of 
Z)W and D^ 2 \ one has a simple, "on-the-fly" criterion to 
quantify the accuracy of this estimation with an empha- 
sis on the precise modeling of entropy-consuming events. 
The method presently used for the verification of D^> 
and Z?( 2 ) involves the numerical solution of the FPE 
with the estimated drift and diffusion coefficients and 
a comparison with the underlying experimental trajecto- 
ries 
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which is, of course, a much more cumbersome 
procedure. 

It is interesting to elucidate some characteristics of 
the entropy-consuming trajectories. To contrast entropy- 
consumption with entropy production, we show in Fig. [3] 
the average of 50 extreme sequences u(r) giving rise to 
very small and very large values of AStot respectively. 
These averages display the distinct features common to 
all individual trajectories of the corresponding class. As 
expected, trajectories giving large and small values of 
AStot look rather different from each other. Large en- 
tropy production, as shown in the bottom panel of Fig. [3J 
is related to a continuous decrease of u for decreasing 
r. Rather similar to macroscopic thermodynamics, dissi- 
pation therefore stems from the transfer of energy from 
large-scale to small-scale motion. In contrast, negative 
values of AStot require violent fluctuations at small scales 
together with a smooth flow at large scales as shown in 
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the top panel of Fig. |3l Therefore, the same class of 
fluctuations that causes small-scale intcrmittency in de- 
veloped turbulence also ensures the good convergence of 
the integral fluctuation theorem ©. 

This connection becomes also apparent when study- 
ing the deviations from dimensional scaling in developed 
turbulence. Consider the moments 



S n {r) = / u n p{u,r)du 



(10) 



of the distribution p(u,r). The self-similar eddy- 
hierarchy suggests scaling laws for these moments of the 



form S n (r) 



defining the scaling exponents Cn- A 



relation for these exponents, the so-called K62 scaling, 
was proposed in 1962 by Kolmogorov and Oboukhov on 
the basis of dimensional analysis and some simplifying as- 
sumptions about the stochastic energy transfer between 



scales 



44, 4 



n n(n — 3) 
C " = 3 ~ M 18 



(11) 



The intermittency factor \i describes deviations from 
pure dimensional (K41 (46|) scaling. It is an experimen- 
tal fit factor with typical values of about 0.25 32| . For 
the data of Fig.[T]/x w 0.227. 
Choosing 



D {1) (u,r) = -- 



D (2) (w,r) 



9r 



A* 2 



(12) 



the stochastic dynamics ((SJ reproduces the K62 scaling 
pip for the moments (jTUJ) for any initial distribution 
P(ul,L) [24|,[38|. Note that this is already the most gen- 
eral case: In order to find a scaling law S n (r) oc r* n 
from the Fokker-Planck dynamics ([2]), one must have 
D<M ~ u/r and ~ u 2 /r 

These dependencies are, however, also special with re- 
spect to the fluctuation theorem (JSJ . Given (fT2j) , we may 
transform to logarithmic "time" log L/r to end up with 
a FPE describing a stationary process without external 
driving. The fluctuation theorem then merely describes 
the relaxation process from an initial non-equilibrium dis- 
tribution to the stationary state p st = 6(u) where all 



S n (r) ->■ [35J, |4j|. Corrections to K62 scaling there- 
fore correspond to a non-trivial "time" dependence of 
drift and diffusion coefficients in the FPE and hence ex- 
press genuine non-equilibrium dynamics along the eddy- 
hierarchy. 

To highlight the sensitivity of the fluctuation theorem 
to small-scale intcrmittency, we specify (j8]) to the drift 
and diffusion coefficients (fT2|) of K62 scaling. Using ((6]) 
and (0, we find 



< Pr(u r ) 

u u L p L {u L ) 



= 1, 



(13) 



with v = ~ 28. This large value of v is consistent 

with the qualitative picture discussed above: Trajecto- 
ries corresponding to large values of AStot have ul > u ri 



whereas those with negative A5 t ot feature ul < u r . 
Using data from numerical simulations of the Langevin 
equation corresponding to (fj"2"j) , we indeed find a smooth 
convergence of (fT3"| for sample sizes of 10 4 or larger. 

The crucial point, however, is that (|13[) fails dramat- 
ically for realistic turbulent flows. Using again the ex- 
perimental data of Fig. [TJ the average in (|13[) results into 
about 10 70 instead of 1! The value 1 is only approached if 
small-scale fluctuations occur with the frequency charac- 
teristic for the K62 model. The much more frequent and 
stronger fluctuations of a realistic turbulent flow, how- 
ever, cause the rapid divergence of (|13p . which we ex- 
plain by the well known fact that K62 underestimates 
the frequency of large fluctuations on small scales (i.e., 
the scaling (fTTj) is only good for n < 10 [H, |49[). Hence, 
the corresponding failure of K62 to accurately describe 
the tails of p(u, r) is most strikingly demonstrated by the 
breakdown of (fl3|) . 

In conclusion, we have shown that the violent small- 
scale fluctuations in turbulent flows make developed tur- 
bulence an interesting model system for stochastic ther- 
modynamics. The integrated fluctuation theorem for en- 
tropy production represents an extremely sensitive in- 
dicator for the precise modeling of small-scale intcrmit- 
tency. Moreover, it may be used as a simple "sum-rule" 
to quantify the accuracy of parameter estimation from 
experimental data drawn from turbulent flows. Also 
other models of developed turbulence l ike those yield- 
ing scaling laws different from K62 [HT 



51|, 



jropagator 

methods [52| , or field-theoretic approaches [4l|, |53j cor- 



respond to a Markovian dynamics of velocity increments 
on the eddy-hierarchy [ID, |4?], [H|. It should be inter- 
esting to apply the presented analysis using fluctuation 
theorems also to these approaches. 
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